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ABSTRACT 

The  Sounding  of  the  Atmosphere  using  Broadband  Emission  Radiometry  (SABER)  experiment  was  launched 
onboard  the  TIMED  satellite  in  December,  2001.  SABER  is  designed  to  provide  measurements  of  temperature, 
constituents,  and  the  key  radiative  and  chemical  sources  and  sinks  of  energy  in  the  mesosphere  and  lower 
thermosphere  (MLT).  SABER  measures  Earth  limb  emission  in  10  broadband  radiometer  channels  ranging 
from  1.27  fim  to  17  //m.  Measurements  are  made  both  day  and  night  over  the  latitude  range  from  52°S  to  83° N 
with  alternating  hemisphere  coverage  every  60  days.  In  this  paper  we  concentrate  on  retrieved  profiles  of  kinetic 
temperature  (Tk)  and  CO2  volume  mixing  ratio  (vmr),  inferred  from  observed  15  //m  and  4.3  fim  limb  emissions. 
SABER-measured  limb  radiances  are  in  non-local  thermodynamic  equilibrium  (non-LTE)  in  the  MLT  region. 
The  complexity  of  non-LTE  radiation  transfer  combined  with  the  large  volume  of  data  measured  by  SABER 
requires  new  retrieval  approaches  and  radiative  transfer  techniques  to  accurately  and  efficiently  retrieve  the  data 
products.  In  this  paper  we  present  the  salient  features  of  the  coupled  non-LTE  Tk/C02  retrieval  algorithm, 
along  with  preliminary  results. 

Keywords:  Remote  sensing,  non-local  thermodynamic  equilibrium  (non-LTE),  thermal  structure,  carbon  diox¬ 
ide,  mesosphere,  lower  thermosphere,  middle  atmosphere 


1.  INTRODUCTION 

On  December  7,  2001,  NASA  launched  the  Sounding  of  the  Atmosphere  using  Broadband  Emission  Radiometry 
(SABER)  experiment  onboard  the  TIMED  satellite.  The  satellite  was  placed  in  a  74.1  degree  inclined,  625 
km  orbit  by  a  Delta  II  rocket.  The  SABER  instrument  is  designed  to  provide  measurements  of  temperature, 
constituents,  and  the  major  radiative  and  chemical  sources  and  sinks  of  energy  in  the  MLT  region.  The  primary 
science  goal  of  SABER  is  to  achieve  major  advances  in  our  understanding  of  the  structure,  energetics,  chemistry, 
and  dynamics  of  the  atmospheric  region  extending  from  60  km  to  180  km  altitude.  This  will  be  accomplished 
using  the  measurement  approach  of  spectral  broadband  emission  radiometry.  SABER  scans  the  horizon  and 
observes  limb  emission  in  10  spectral  channels  ranging  from  1.27  fim  to  17  [im.  The  observed  limb  emission 
profiles  are  analyzed  to  provide  vertical  profiles,  with  approximately  2  km  altitude  resolution,  of  the  following 
parameters:  Tk;  CO2,  O3  and  H2O  vmr;  atomic  oxygen  and  atomic  hydrogen;  volume  emission  rates  due  to 
O2OA),  OH(i/=3,4,5),  OH(i/=7,8,9),  and  NO  at  5.3  jum;  key  atmospheric  cooling  rates,  solar  heating  rates, 

*  Correspondence:  E-mail:  c.j.mertens@larc.nasa.gov;  Telephone:  757-864-2179;  Fax:  757-864-6326 
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chemical  heating  rates,  and  geostrophic  winds.  Measurements  are  made  both  day  and  night  over  the  latitude 
range  from  53°S  to  82°N  with  alternating  hemisphere  coverage  every  60  days. 

In  the  MLT  region,  SABER  measures  in  each  channel  -  with  the  exception  of  the  1.27  fim  channel,  which 
mainly  measures  emission  from  the  C^1  A)  electronic  transition  -  infrared  radiation  emitted  by  molecules  whose 
vibration-rotation  bands  are  in  non-local  thermodynamic  equilibrium  (non-LTE).  Non-LTE  kinetic  processes 
and  non-LTE  radiation  transfer  complicate  the  analysis  and  retrieval  of  the  SABER  data  products,  as  well  as 
impose  a  significant  computational  burden  on  data  processing.  New  retrieval  approaches  and  radiation  transfer 
techniques  are  required  to  accurately  and  efficiently  retrieve  the  data  products  from  the  large  volume  of  SABER 
non-LTE  emission  measurements.  In  particular,  the  focus  of  this  paper  is  on  the  simultaneous  non-LTE  kinetic 
temperature  (Tk)  and  CO2  retrieval  approach. 

Kinetic  temperature  is  retrieved  from  SABER’s  broadband  measurement  of  CO2  15  /im  limb  emission.  This 
technique  was  developed  more  than  30  years  ago.1  In  these  early  experiments,  a  basic  assumption  was  that  CO2 
was  well  mixed  and  its  volume  mixing  ratio  well  known.  Another  key  assumption  was  that  the  observed  CO2 
vibration-rotation  bands  were  in  LTE.  These  assumptions  were  sufficient  for  previous  sensors  whose  sensitivity 
did  not  permit  limb  radiance  measurements  much  above  70  km  tangent  height. 

The  assumption  of  LTE  in  the  CO2  infrared  bands  and  the  assumption  of  uniformly  mixed  and  well  known 
CO2  vmr,  described  in  the  last  paragraph,  are  not  valid  in  the  MLT  region.  The  variability  of  CO2  vmr  requires 
that  Tk  and  CO2  vmr  be  retrieved  simultaneously.  SABER  measures  CO2  limb  emission  in  the  15  /im  spectral 
interval  to  approximately  120  km  in  altitude  for  the  purpose  of  determining  Tk.  SABER  also  observes  4.3  /im 
CO2  limb  emission  to  over  160  km  altitude  during  the  day  and  to  approximately  130  km  at  night.  Measurements 
of  CO2  15  /im  and  CO2  4.3  //m  limb  radiance  are  combined  to  retrieve  Tk  and  CO2  vmr  in  the  MLT  region. 
Non-LTE  processes  are  rigorously  accounted  for  in  the  retrieval  scheme. 

Derived  profiles  of  Tk  in  the  MLT  are  necessary  to  understand  the  thermal  structure  of  this  region,  one  of 
SABER’s  primary  goals.  Accurate  knowledge  of  Tk  and  CO2  are  necessary  to  quantify  the  radiative  cooling  of 
this  region  of  the  atmosphere,  largely  dominated  by  infrared  emission  from  the  CO2  15  /im  bands.  Moreover,  Tk 
is  a  key  input  into  the  retrieval  of  SABER’s  other  data  products.  The  global  distribution  of  CO2  is  interesting 
in  its  own  right,  as  its  latitudinal  and  seasonal  variations  are  not  well  known  and  observed,  and  the  mechanism 
for  its  departure  from  its  uniformly  mixed  value  is  not  well  characterized. 

In  a  previous  work,  Mertens  et  al.2  presented  a  non-LTE  Tk  retrieval  algorithm  based  on  observations  of  CO2 
15  /im  broadband  lirfib  emission  measurements,  assuming  CO2  abundance  was  known.  They  demonstrated  the 
algorithm  using  model  atmospheres,  and  studied  the  sensitivity  of  retrieved  Tk  to  key  atmospheric  and  kinetic 
parameters  used  in  the  non-LTE  CO2  model.  One  important  conclusion  from  this  work  was  that  CO2  vmr 
needed  to  be  known  to  15%  or  better  in  order  to  achieve  SABER’s  goal  of  retrieving  Tk  to  3  K  or  better 
below  100  km.  As  a  result  of  this  study,  the  non-LTE  Tk  retrieval  algorithm  of  Mertens  et  al.  has  been 
expanded  and  updated  to  enable  a  simultaneous  non-LTE  retrieval  of  Tk  and  CO2.  In  this  paper,  we  present 
the  salient  features  of  the  coupled  non-LTE  Tk/C02  retrieval  algorithm.  The  algorithm  is  demonstrated  by 
showing  preliminary  retrievals  from  SABER  measurements. 

2.  RETRIEVAL  APPROACH 

Kinetic  temperature  is  retrieved  in  the  stratosphere  using  SABER  measured  radiance  from  two  CO2  15  /im 
channels,  a  narrow  bandpass  channel  (650-695  cm-1)  and  a  wide  bandpass  channel  (580-760  cm-1).  The 
two  CO2  channels  are  used  to  register  pressure  with  altitude  in  the  stratosphere  and  infer  Tk  assuming  LTE 
conditions  and  assuming  CO2  is  uniformly  mixed  and  known.  This  approach  is  similar  to  the  two-color  technique 
described  in  Ref..1  The  LTE  assumption  breaks  down  in  the  mesosphere  for  the  infrared  CO2  bands,  and  CO2 
vmr  is  no  longer  uniformly  mixed.  The  simultaneous  non-LTE  Tk/CC>2  retrieval  algorithm  is  then  employed  to 
infer  Tk  and  CO2  vmr  in  the  MLT  using  measured  radiance  from  the  CO2  15  //m  narrow  channel  and  the  CO2 
4.3  //m  channel  (2320-2400  cm-1). 

The  LTE-retrieved  Tk  and  pressure  described  in  the  preceding  paragraph  provide  the  lower  boundary  con¬ 
ditions  for  the  non-LTE  Tk/CC>2  retrieval.  The  lower  boundary  condition  in  the  CO2  vibrational  temperature 
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model  (described  in  the  next  section)  is  that  the  source  function  in  the  equation  of  radiative  transfer  be  given 
by  the  Planck  function,  requiring  that  the  C02  bands  be  optically  thick  and  in  LTE.  In  the  retrieval  inversion 
approach,  pressure  (p)  is  obtained  by  vertically  integrating  the  barometric  equation  from  the  lower  boundary 
altitude,  requiring  that  the  LTE-retrieved  Tk(p)  be  accurate  at  the  lower  boundary  altitude.  Taking  these 
factors  into  consideration,  the  lower  boundary  altitude  of  the  non-LTE  Tk/C02  retrieval  is  nominally  taken  to 
be  40  km. 

The  non-LTE  Tk/C02  retrieval  model  is  comprised  of  two  main  components:  (1)  the  forward  radiance  model 
and  (2)  the  inversion  model.  These  two  components  are  described  in  more  detail  below. 

2.1.  Forward  Model 

The  forward  radiance  model  of  the  retrieval  algorithm  is  the  component  that  simulates  the  measured  radiance 
along  the  limb  line-of-sight.  The  forward  model  itself  is  composed  of  two  parts:  (1)  the  vibrational  temperature 
model  (Tv)  and  (2)  the  limb  radiance  model.  Limb  radiance  is  calculated  using  BANDPAK,3  now  expanded 
for  application  to  non-LTE  calculations.4  The  non-LTE  formulation  in  BANDPAK  is  a  broadband  extension  of 
the  line-by-line  approach  described  by  Edwards  et  al5  and  initially  demonstrated  by  Mlynczak  et  al.6  There  are 
eighteen  C02  15  /im  bands  and  one  O3  14.1  /zm  band  that  contribute  to  the  total  limb  radiance  in  the  SABER 
C02  15  /2m  narrow  channel  spectral  bandpass,  above  the  lower  boundary  of  the  non-LTE  Tk/C02  retrieval 
model.  The  SABER  4.3  //m  channel  has  seventeen  C02  bands  that  contribute  to  the  total  limb  radiance  at 
altitudes  of  approximately  70  km  and  above.  Vibrational  temperatures  for  the  nineteen  bands  that  emit  in  the 
C02  15  /xm  narrow  channel  bandpass  and  the  seventeen  bands  that  emit  in  the  C02  4.3  /zm  channel  bandpass 
comprise  the  non-LTE  inputs  into  the  limb  radiance  model. 

The  nineteen  bands  that  contribute  to  the  total  limb  radiance  in  the  SABER  C02  15  /im  narrow  channel 
bandpass  can  be  grouped  into  seven  band-groups:  the  fundamental  v2  band  of  the  major  (626)  isotope  (010- 
626);  the  fundamental  v2  bands  of  the  minor  (636,  628,  and  627)  isotopes  (010-MIN);  the  first  i/2  hot  bands  of 
the  626  isotope  ({020}-626);  the  first  v2  hot  bands  of  the  minor  isotopes  ({020}-MIN);  the  second  i/2  hot  bands 
of  the  626  isotope  ({030}-626);  the  O3  major  isotopic  v2  fundamental  band;  and  the  ‘remaining’  C02  bands 
that  contribute  to  the  C02  15  /im  narrow  channel  limb  radiance.  In  the  above  notation,  (0i/20)  refers  to  a  single 
transition  and  {0j/20}  refers  to  a  group  of  transitions.  Panels  (a)  and  (b)  of  Figure  1  show  SABER  C02  15 
/im  narrow  channel  limb  radiance  simulations  for  the  US  Standard  model  atmosphere:  limb  radiance  from  the 
seven  band-groups  listed  above,  the  total  limb  radiance,  the  contribution  of  each  band-group  to  the  total  limb 
radiance,  and  the  noise  level  (NER:  noise  equivalent  radiance)  in  the  SABER  C02  15  /z m  narrow  channel.  The 
C02  ‘remainder’  band  contains  contributions  from  all  other  C02  bands  found  on  the  HITRAN  1996  database 
that  are  not  specified  in  the  band-groups  above.  Although  there  are  no  significant  non-LTE  effects  in  the  limb 
radiance  from  these  ‘remainder’  bands,  their  cumulative  contribution  to  the  total  limb  radiance  is  significant  and 
must  be  modeled,  as  Figure  1  indicates.  The  C02  ‘remainder’  band  is  treated  as  a  pseudo- vibrational  band  -  the 
non-LTE-to-LTE  vibrational  state  population  ratios  (see  Edwards  et  al.5)  are  set  equal  to  the  non-LTE-to-LTE 
population  ratios  of  the  626  first  i/2  hot  band.  This  step  eliminates  anomalous  thermospheric  contributions  and 
guarantees  that  the  correct  LTE  contributions  are  included  at  altitude  below  approximately  70  km,  where  their 
‘true’  contribution  becomes  important,  as  Figure  1  indicates. 

The  bands  that  contribute  most  to  the  SABER  C02  15  /zm  narrow  channel  spectral  bandpass  are  the  010- 
626,  {020}-626,  and  010-MIN  band-groups.  The  010-626  band  dominates  the  limb  radiance  except  near  80  km, 
where  the  {020}-626  and  010-MIN  bands  rival,  if  not  exceed,  the  contribution  due  to  010-626.  The  contribution 
from  the  other  band-groups  are  much  smaller  than  those  mentioned  above;  however,  they  must  be  rigorously 
modeled  in  order  to  retrieve  Tk  accurately. 

Analogous  to  the  discussion  of  the  bands  that  emit  in  the  SABER  C02  15  /zm  narrow  channel  bandpass, 
the  seventeen  bands  that  contribute  to  the  total  limb  radiance  in  the  C02  4.3  /zm  channel  can  be  grouped  into 
six  band-groups:  the  fundamental  v 3  band  of  the  626  isotope  (001-626);  the  fundamental  2/3  bands  of  the  minor 
isotopes  (001-MIN);  the  first  4.3  /zm  hot  band  of  the  626  isotope  (011-626);  and  the  second,  third  and  fourth 
4.3  /zm  hot  bands  of  the  626  isotope  ({021}-626,  {031}-626,  and  {041}-626)).  Panels  (c)  and  (d)  of  Figure  1 
for  the  SABER  C02  4.3  /zm  channel  are  analogous  to  panels  (a)  and  (b)  for  the  C02  15  /zm  narrow  channel. 
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Figure  1.  Simulations  of  daytime  (solar  zenith  angle  (SZA)  60  degrees)  non-LTE  line-of-sight  (LOS)  limb  radiance  for 
the  SABER  CO2  15  /zm  narrow  channel  (panels  (a)  and  (b))  and  CO2  4.3  /*m  channel  (panels  (c)  and  (d))  for  the  US 
Standard  atmosphere.  See  text  for  details. 

The  bands  that  contribute  most  to  the  CO2  4.3  /zm  channel  radiance  are  the  001-626,  001-MIN,  011-626,  and 
{021  }-626  band-groups.  Above  110  km,  limb  emission  in  the  4.3  /zm  channel  is  dominated  by  the  001-626  band. 
Between  95  km  and  110  km,  the  four  major  band-groups  mentioned  above  are  comparable  in  their  contribution 
to  the  total  limb  radiance.  Below  95  km,  limb  emission  is  dominated  by  the  {021  }-626  band-group.  However, 
similar  to  the  case  for  the  SABER  CO2  15  /zm  narrow  channel,  all  band-groups  shown  in  panels  (c)  and  (d) 
must  be  rigorously  modeled  in  order  to  accurately  retrieve  CO2  vmr. 

Figure  2  shows  a  model  calculation  of  the  CO2  Tv’s  for  the  US  Standard  atmosphere.  The  solar  pumped 
states  are  shown  in  panels  (a)-(d).  The  15  /zm  Tv’s  are  shown  in  panel  (e).  The  solar  pumped  states  are 
responsible  for  the  emission  in  the  SABER  CO2  4.3  /zm  channel.  Furthermore,  they  have  an  indirect  effect  on 
emission  in  the  CO2  15  /zm  channel  through  collisional  and  radiative  coupling  to  the  CO2  ^2  vibrational  state 
manifold,  as  indicated  in  the  15  /zm  day /night  differences  shown  in  panel  (f).  The  vibrational  temperatures 
are  calculated  from  the  operational  CO2  Tv  model,  which  is  based  on  the  Modified  Curtis  Matrix  approach 
advanced  by  Lopez-Puertas  et  al..7  The  operational  CO2  Tv  model  uses  BANDPAK  to  perform  all  the  radiation 
transfer  calculations  and  will  be  described  in  more  detail  in  Refs.4  and.8 

2.2.  Inversion  Method 

Kinetic  temperature  and  CO2  vmr  are  retrieved  by  successively  iterating  between  two  independent  retrieval 
modules:  one  for  Tk(p),  assuming  CO2  is  known,  and  the  other  one  for  CO2  vmr,  assuming  Tk(p)  is  known. 
The  iteration  stops  when  successive  CO2  profiles  have  relaxed  sufficiently  over  the  altitude  range  of  the  non-LTE 
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Figure  2.  Simulations  of  daytime  (solar  zenith  angle  (SZA)  60  degrees)  C02  vibrational  temperatures  for  US  Standard 
atmosphere.  The  following  states  are  shown:  (a)  {041},  (b)  {031},  (c)  {021},  (d)  Oil  and  001,  (e)  {0i/20},  (f)  {0i/20} 
day /night  differences.  Individual  states  are  labeled  in  HITRAN  database  notation. 


retrieval  model.  The  global  convergence  criteria  for  the  coupling  between  the  Tk  and  C02  retrievals  is  that 
the  average  difference  between  two  successive  C02  vmr  retrievals,  i.e.,  averaged  over  all  the  retrieved  tangent 
levels,  does  not  exceed  a  user-specified  percentage.  Retrieval  simulations  for  various  model  atmospheres  (US 
Standard,  polar  summer,  and  polar  winter)  have  shown  that  the  differences  between  the  last  and  next-to-the-last 
retrieved  C02  profiles  are  small  enough  not  to  significantly  affect  an  additional  Tk  retrieval;  thus,  justifying  the 
termination  of  the  non-LTE  Tk/C02  retrieval  algorithm  based  on  the  convergence  of  the  C02  profile.  Below 
we  briefly  describe  the  separate  Tk  and  C02  retrieval  modules.  A  detailed  algorithm  description  will  appear 
elsewhere.9 

First  we  describe  the  non-LTE  Tk  retrieval  module.  There  are  two  primary  relaxation  loops.  In  the  inner 
loop,  a  Tk  profile  is  retrieved  using  the  onion-peel  approach10  while  pressure,  C02  vmr,  and  the  Tv’s  are  held 
fixed.  Kinetic  temperature  is  retrieved  at  each  tangent  altitude  by  adjusting  the  local  Tk  until  the  modeled 
radiance  matches  the  measured  radiance  within  the  convergence  criterion.  The  temperature  is  adjusted  using 
the  Levenberg-Marquardt  method.11  The  inner  loop  convergence  criterion  is  a  requirement  that  the  difference 
between  successive  iterations  is  much  less  than  the  expected  errors  in  the  solution.12  The  fraction  of  the 
solution  error  required  to  satisfy  the  convergence  criterion  is  user-specified. 

The  onion-peel  approach  is  critical  to  retrievals  from  C02  limb  emission  measurements  in  the  MLT  region, 
since  the  C02  15  fim  limb  radiance  for  mesospheric  tangent  heights  is  dominated  by  emission  from  higher 
altitude  layers.13  The  same  is  also  true  for  C02  4.3  /im  limb  radiance.  The  onion-peel  technique  ensures  that 
the  modeled  emission  matches  the  measured  radiance  from  the  upper  altitude  layers,  even  though  the  retrieved 
temperature-pressure-C02  combination  may  be  incorrect  at  intermediate  steps  in  the  relaxation  process.  For  a 
particular  limb  path,  the  effect  is  greater  sensitivity  to  the  local  Tk  at  the  sought-after  tangent  altitude. 
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In  the  outer  relaxation  loop,  the  pressure  profile  is  rebuilt  from  the  lower  boundary  altitude  using  the  onion- 
peel  retrieved  Tk  profile  and  the  barometric  pressure  law.  The  vibrational  temperatures  are  updated  using 
the  CO2  Tv  model  with  the  previously  retrieved  Tk,  CO2  vmr,  and  pressure  profile  as  input.  The  onion-peel 
retrieval  (inner  loop)  is  repeated  until  the  entire  inferred  Tk  profile  relaxes  within  the  convergence  criterion,  a 
criterion  similar  to  the  inner  loop  convergence  criterion  described  above. 

Next  we  describe  the  non-LTE  CO2  vmr  retrieval  module.  Because  of  the  severe  nonlinearities  in  the 
radiation  transfer  along  the  limb  line-of-sight  in  the  SABER  CO2  4.3  fim  channel  bandpass,  the  CO2  vmr 
retrieval  module  is  composed  of  a  juxtaposition  of  two  retrieval  approaches.  The  first  approach  is  the  nonlinear 
relaxation  method  of  Twomey-Chahine,14  modified  here  for  limb  path  geometry  and  broadband  radiance 
measurements.  In  the  current  approach,  the  kernel  function,  or  weighting  function,  is  the  contribution  of  each 
layer  along  the  limb  line-of-sight  to  the  total  limb  radiance,  divided  by  the  layer-averaged  CO2  abundance. 
Similar  to  the  non-LTE  Tk  retrieval,  the  CO2  vmr  is  retrieved  at  each  tangent  altitude  using  the  onion-peel 
method.  The  convergence  criterion  was  described  in  the  first  paragraph  of  this  section.  If  the  criterion  is  not 
satisfied,  the  non-LTE  Tk  retrieval  is  repeated,  followed  by  a  subsequent  CO2  vmr  retrieval  using  the  modified 
Twomey-Chahine  approach.  Once  the  global  convergence  criterion  is  satisfied,  one  final  CO2  vmr  retrieval 
is  performed  using  the  Levenberg-Marquardt  approach,  analogous  to  the  inner  loop  relaxation  scheme  of  the 
non-LTE  Tk  retrieval  described  above. 

The  nonlinearities  in  the  CO2  4.3  /zm  radiation  transfer  are  too  severe  to  use  a  retrieval  algorithm  based 
on  Newtonian  iteration,  applicable  to  weakly  nonlinear  problems,15  or  a  Levenberg-Marquardt  approach, 
applicable  to  moderately  nonlinear  problems.15  Consequently,  the  nonlinear  relaxation  method  of  Twomey- 
Chahine  is  used  to  infer  CO2  vmr  in  the  coupled  non-LTE  Tk/CC>2  retrieval  algorithm.  However,  optimal 
estimation  approaches,  such  as  the  Levenberg-Marquardt  approach,  provide  useful  statistical  information:  the 
estimated  solution  error,  quality  metrics  such  as  the  x2,  or  penalty  function,  and  more  theory-based  convergence 
criterion.  Therefore,  the  final  CO2  vmr  retrieval  uses  the  Levenberg-Marquardt  algorithm  to  refine  the  modified 
Twomey-Chahine  approach  and  provide  the  useful  statistical  information  described  above. 

3.  RESULTS  AND  DISCUSSION 

In  this  section  we  present  preliminary  non-LTE  Tk/CC>2  retrievals  from  SABER  daytime  observations.  At  the 
present  time,  CO2  is  not  being  retrieved  at  night.  After  launch  it  was  realized  that  the  noise  level  in  the  4.3 
fim  channel  was  much  lower  than  anticipated,  offering  the  potential  to  retrieve  CO2  at  night.  We  are  currently 
investigating  whether  there  is  enough  information  content  in  the  4.3  /zm  measurements  in  the  absence  of  solar 
radiative  excitation  to  retrieve  CO2  at  night.  In  the  short  term,  during  nighttime  observations,  we  present 
non-LTE  Tk  retrievals,  with  the  CO2  profile  taken  from  the  TIME-GCM  climatology.  There  are  a  number  of 
additional  atmospheric  parameters  required  as  input  into  the  CO2  Tv  model:  O2,  N2,  0(3P),  and  0(1D),  for 
example.  In  the  preliminary  retrieved  profiles  shown  in  this  section,  the  atmospheric  parameters  listed  above 
were  obtained  from  the  TIME-GCM  climatology  produced  for  SABER  analysis. 

Figure  3  shows  nighttime  non-LTE  Tk  retrievals  on  March  3,  2002,  near  69°N.  The  retrieved  Tk  profiles  are 
compared  to  the  Liibken  and  von  Zahn16  climatological  monthly  mean  Tk  profile  for  March  at  69°N.  The  monthly 
mean  profile  was  generated  from  Tk  measurements  taken  from  sodium  Lidar,  falling  spheres,  and  rocketborne 
mass  spectrometer  and  ionization  gauge  measurements  (see  Ref.16  for  details).  Between  60  km  and  105  km, 
the  differences  between  the  individual  SABER-derived  Tk  profiles  and  the  monthly  mean  profile  are  within  the 
natural  variability  in  the  monthly  mean  Tk  profile  for  March  at  69°N.16  Below  60  km  the  SABER  Tk  profiles 
are  substantially  cooler  than  the  climatological  mean;  above  105  km  the  SABER  Tk  profiles  are  substantially 
warmer  than  the  climatological  mean.  Further  analysis  is  required  to  understand  these  differences. 

The  following  figures  demonstrate  the  coupled  non-LTE  Tk/CC>2  retrieval  algorithm  for  daytime  SABER 
observations.  Figure  4  shows  retrieved  CO2  profiles  on  January  28  and  March  3,  2002.  For  each  day,  retrieved 
profiles  are  shown  at  latitudes  of  approximately  44° N  and  27° S.  The  retrieved  profiles  are  shown  as  solid 
lines.  Climatological  mean  TIME-GCM  model17  profiles  which  most  closely  correspond  to  the  time,  season  and 
geo-location  of  the  SABER  measurements  are  shown  as  dashed  lines. 
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Figure  3.  Nighttime  non-LTE  Tk  retrievals  from  SABER  measurements  on  March  3,  2002,  near  69°N.  The  latitude, 
longitude,  orbit  and  scan  numbers,  respectively,  for  the  SABER  measurements  are  indicated  in  the  legend.  The  retrieved 
Tk  profiles  are  compared  to  the  climatological  monthly  mean  March  Tk  profile  from  Liibken  and  von  Zahn.16 


For  equinox  conditions,  model  simulations  indicate  that  the  latitudinal  distribution  of  CO2  is  symmetrically 
distributed  about  the  equator  (see  Ref.18  and  references  therein).  Consequently,  the  CO2  profiles  for  March 
3  at  44°N  and  27°S  should  be  very  similar,  as  Figure  4  indicates.  However,  in  some  cases,  SABER-derived 
CO2  profiles  suggest  a  much  larger  depletion  of  CO2  in  the  upper  mesosphere  and  lower  thermosphere  than 
predicted  by  the  TIME-GCM  model  and  previous  rocket  borne  mass  spectrometer  measurements,18  similar  to 
the  findings  of  other  recent  CO2  retrievals  from  infrared  emission  measurements  (see  Refs.19-21).  From  Figure  4, 
the  TIME-GCM  simulations  predict  that  for  equinox  conditions,  the  CO2  abundance  for  27° S  is  greater  than 
the  CO2  abundance  for  44° N  above  100  km.  The  SABER  retrievals  show  just  the  opposite  effect:  the  CO2 
abundance  at  27°S  is  significantly  depleted  relative  to  the  CO2  abundance  at  44°N.  However,  these  differences 
could  be  due  to  the  fact  that  we  are  comparing  TIME-GCM  equinox  profiles  with  SABER  observations  on  one 
day  in  March,  which  may  not  be  representative  of  equinox  conditions. 

Model  simulations  indicate  that,  for  solstice  conditions,  CO2  vmr  remains  constant  up  to  higher  altitudes  in 
the  summer  hemisphere  and  departs  from  its  constant  value  at  lower  altitudes  in  the  winter  hemisphere,  an  effect 
governed  by  upward  transport  in  the  summer  hemisphere  and  downward  transport  in  the  winter  hemisphere.18 
SABER  CO2  retrievals  support  this  conclusion,  as  indicated  in  Figure  4  for  the  January  28  profiles.  However, 
the  retrieved  CO2  profile  at  44°N  is  significantly  depleted  with  respect  to  the  corresponding  TIME-GCM  profile, 
while  the  corresponding  TIME-GCM  profile  for  the  27°S  measurement  seems  to  fit  the  measurement  scenario 
quite  well. 

Figure  5  shows  the  Tk  profiles  that  were  retrieved  simultaneously  with  the  CO2  profiles  shown  in  the  previous 
figure.  Similar  to  Figure  4,  retrieved  profiles  are  shown  on  January  28  and  March  3,  2002,  at  approximately 
44°N  and  27°S.  Retrieved  profiles  are  shown  as  solid  lines  and  TIME-GCM  climatological  mean  Tk  profiles  are 
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Figure  4.  Daytime  non-LTE  Tk/C02  retrievals  from  SABER  measurements  on  January  28  and  March  3,  2002.  In  this 
figure  we  show  the  retrieved  CO2  profiles.  The  month,  day,  latitude,  longitude,  orbit  and  scan  numbers,  respectively, 
for  the  SABER  measurements  are  indicated  in  the  legend.  The  retrieved  CO2  profiles  are  compared  to  the  TIME-GCM 
climatological  mean  that  most  closely  corresponds  to  the  SABER  measurement  time,  season  and  geo-location.  The 
TIME-GCM  profiles  shown  were  computed  for  two  time  periods:  daytime  winter  solstice  (indicated  by  SOL.D355)  and 
daytime  spring  equinox  (indicated  by  EQNJD080). 


SABER-Derived  Kinetic  Temperature  Profiles  for  JAN28/MAR03  2002 


Figure  5.  Daytime  non-LTE  Tk/CC>2  retrievals  from  SABER  measurements  on  January  28  and  March  3,  2002.  See 
caption  to  Figure  4  for  details. 
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shown  as  dashed  lines. 

The  most  noticeable  features  in  Figure  5  above  110  km  are  found  in  the  Tk  profiles  for  January  28  and 
March  3  at  44°N.  There’s  an  inversion  layer  in  the  March  3  profile  and  Tk  decreases  with  increasing  altitude  in 
the  January  28  profile.  Neither  of  these  features  are  real;  rather,  they’re  due  to  large  noise  spikes  in  the  CO2 
15  /im  narrow  channel  radiance  profiles.  In  general,  measurement  noise  will  prohibit  reasonable  single  profile 
Tk  retrievals  much  above  110  to  115  km.  However,  daily  and  weekly  averaging  of  radiance  profiles  before  the 
retrievals  should  significantly  extend  the  altitude  range. 

The  SABER-derived  Tk  profiles  in  Figure  5  show  a  great  deal  of  vertical  wave  structure,  as  do  the  SABER 
profiles  in  Figure  3,  suggestive  of  rather  strong  tidal  signatures.  The  SABER  Tk  profiles  in  Figure  5  are  generally 
warmer  than  the  TIME-GCM  profiles  between  70  km  and  90  km.  In  particular,  the  January  28  profile  at  27°S 
and  the  March  3  profile  at  44°N  are  significantly  warmer  than  the  model  profiles  between  70  km  and  90  km. 
On  the  other  hand,  the  SABER  Tk  profile  for  March  3  at  27°S  is  reasonably  close  to  the  corresponding  TIME- 
GCM  profile  below  110  km.  The  SABER  January  28  profile  at  44°N  has  a  double  temperature  minimum  at 
approximately  93  km  and  110  km.  The  corresponding  model  profile  has  a  double  temperature  minimum  at 
approximately  75  km  and  100  km.  For  all  the  profiles  shown  in  Figure  5,  the  SABER  Tk  profiles  are  cooler 
than  the  TIME-GCM  profiles  below  approximately  65  km. 

4.  CONCLUSIONS 

In  this  paper  we  briefly  described  the  coupled  non-LTE  Tk/C02  retrieval  algorithm  used  to  infer  Tk  and 
C02  vmr  from  limb  emission  measurements  observed  from  SABER’s  C02  15  fim  and  C02  4.3  /zm  broadband 
radiometer  channels.  We  have  shown  preliminary  Tk/C02  profiles  from  SABER  observations  that,  at  this  stage 
of  data  processing,  look  reasonable.  Nighttime  Tk  retrievals  near  69°  N  generally  compare  favorably  with  the 
Liibken  and  von  Zahn  climatology.  Daytime  Tk  retrievals  compare  reasonably  well  with  TIME-GCM  profiles, 
although  some  of  the  SABER  Tk  profiles  seem  to  be  modulated  by  quite  large  tidal  signatures.  Retrieved  CO2 
vmr  profiles  seem  to  confirm  that  models  (e.g.,  TIME-GCM)  tend  to  over-predict  C02  in  the  MLT  region. 

At  the  time  of  this  writing,  a  number  of  instrument  corrections  are  being  made  to  improve  the  overall 
quality  of  the  SABER  radiance  measurements.  The  non-LTE  Tk/C02  retrieval  algorithm  presented  here  is 
currently  being  interfaced  and  ingested  into  the  overall  SABER  software  system.  In  the  near  future,  a  number 
of  important  atmospheric  input  parameters  (i.e.,  input  into  the  C02  Tv  model)  will  be  derived  below  100 
km  from  SABER  observations  -  for  example,  0(3P)  and  0(1D).  Thus,  the  excellent  quality  of  the  SABER 
measurements  combined  with  the  non-LTE  Tk/C02  retrieval  algorithm  presented  here,  along  with  the  other 
SABER  non-LTE  retrieval  algorithms,  offer  the  potential  to  significantly  improve  our  understanding  of  the  MLT 
thermal  structure,  chemistry,  and  energetics. 
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